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INTRODUCTION 

The  state  of  the  atmosphere  affects  tactical  Army  operations  and  the 
contributing  weapons  systems  employment.  To  increase  the  friendly 
forces'  relative  combat  power,  meteorological  parameters  must  be  meas- 
ured at  the  time  and  location  of  the  action.  An  important  parameter 
in  atmospheric  measurements,  particularly  for  ballistic  weapon  employ- 
ment, is  the  projectile's  trajectory  integrated  path  crosswinds. 

The  purpose  of  this  report  is  to  present  the  results  of  the  evaluation 
of  the  Saturation  Resistant  Crosswind  Sensor  (SRCS),  a sensor  designed 
to  measure  average  path  crosswinds.  This  evaluation. is  based  on  compara- 
tive data  taken  during  the  test  period  of  31  January  to  8 March  1978  at 
Biggs  Optical  Range  (BOR),  Biggs  Army  Airfield,  Fort  Bliss,  Texas. 

Included  as  part  of  the  report  are  daily  weather  summaries  of  atmospheric 
parameters  prevailing  at  BOR. 

The  SRCS  is  an  active  optical  system  designed  to  measure  the  average 
crosswind  along  a path  from  the  receiver  to  the  transmitter.  Although 
the  SRCS  was  tested  under  actual  field  conditions,  the  evaluation  tests 
were  conducted  cognizant  that  the  SRCS  is  presently  a research  instru- 
ment and  not  yet  intended  for  prolonged  field  use  without  proper  recon- 
figuration. 

This  report  presents  results  of  collected  data  and  an  evaluation  and 
analysis  that  determine  the  accuracy,  reliability,  and  applicability  of 
the  SRCS. 

This  work  was  accomplished  by  personnel  of  the  Atmospheric  Sensing  Divi- 
sion, Atmospheric  Sciences  Laboratory  under  DA  Task  1L162111AH71A3. 

INSTRUMENTATION  REQUIREMENT 

Crosswinds  along  a ballistic  projectile  trajectory  contribute  signifi- 
cantly to  the  total  weapon  error.  Walters^  has  shown  that  direct  fire 
crosswind  errors  on  representative  armor  projectiles  are  significantly 
greater  than  head  and  tail  wind  errors.  To  increase  the  first-round-hit 
probability,  crosswinds  must  be  accurately  known  immediately  before 
firing.  Knowledge  and  application  of  crosswind  information  to  fire  con- 
trol systems  can  also  increase  the  standoff  range  of  friendly  weapons 
without  degrading  weapon  accuracy. 

Several  remote  crosswind  sensors  have  been  developed  in  the  recent  past. 
Four  systems  were  evaluated  at  BOR  during  the  test  period.  The  evaluation 


^D.  L.  Walters,  1975,  "Crosswind  Weighting  Functions  for  Direct-Fire 
Projectiles,"  ECOM  Report  5570,  Atmospheric  Sciences  Laboratory,  White 
Sands  Missile  Range,  NM 
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results  of  the  SRCS  are  presented  in  this  report.  The  results  of  the 
other  system  evaluations  are  reported  separately.^ 

The  SRCS,  an  active  bistatic  sensor  that  measures  weighted  average  cross- 
winds,^  can:  (1)  be  used  in  research  that  may  improve  present  wind  meas- 
urement techniques,  (2)  aid  in  the  characterization  of  the  atmosphere, 
especially  as  related  to  high  energy  laser  propagation  effects,  and  (3) 
serve  as  a baseline  system  in  evaluating  and  analyzing  future  wind  meas- 
urement systems. 

A previous  experimental  prototype  of  the  SRCS‘*»^  proved  the  concept  of 
operational  feasibility,  but  was  unable  to  operate  under  conditions  of 
high  integrated  refractive  index  turbulence  which  can  occur  on  clear 
sunny  days  with  path  lengths  over  500  m and  less  than  10  m above  ground. 
In  1977,  an  exploratory  development  prototype  of  the  improved  sensor  was 
completed.  This  is  the  system  evaluated  and  discussed  here.  The  evalua- 
tion of  this  sensor  will  contribute  to  the  necessary  data  base  for  con- 
tinued future  crosswind  sensor  development  to  satisfy  stated  tactical 
requirements. 


SYSTEM  DESCRIPTION 

The  SRCS  consists  of  a transmitter  and  a receiver  placed  at  opposite  ends 
of  the  path  to  be  measured.  The  transmitter  consists  of  a 55  W quartz- 
iodine  lamp  located  at  the  focus  of  a Fresnel  lens,  and  a power  supply 
for  the  lamp.  The  receiver  houses  an  optical  section  (consisting  of  a 
pair  of  photodiode  detectors  and  situated  close  to  the  focal  points  of 
two  13  cm  Fresnel  lenses)  and  an  electronic  section  that  conditions  the 
signals  from  the  photodiodes  and  determines  the  slope  of  their  covariance. 
Table  1 summarizes  the  SRCS  characteristics,  and  figures  1,  2,  and  3 show 
the  unit  that  was  evaluated. 


2R.  Rodriguez,  1979,  "Evaluation  of  the  Passive  Remote  Crosswind  Sensor," 
ERADCOM  Report,  ASL-TR-0032,  US  Army  Atmospheric  Sciences  Laboratory, 
White  Sands  Missile  Range,  NM 

^G.  R.  Ochs,  6.  F.  Miller,  and  E.  J.  Goldenstein,  1976,  "A  Saturation- 
Resistant  Optical  System  for  Measuring  Average  Wind,"  ECOM  76-2,  Atmo- 
spheric Sciences  Laboratory,  White  Sands  Missile  Range,  NM 

“^R.  S.  Lawrence,  G.  R.  Ochs,  and  S.  F.  Clifford,  1972,  "The  use  of  scin- 
tillations to  measure  average  wind  across  a light  beam,"  Appl  Opt  11, 

No.  2,  239-243 

^G.  R.  Ochs  and  G.  F.  Miller,  1973,  "The  NOAA  optical  system  for  meas- 
uring average  wind,"  NOAA  Technical  Memorandum  ERL  WPL-9 
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TABLE  1.  SATURATION  RESISTANT  CROSSWIND  SENSOR  CHARACTERISTICS 


Receiver 


Function 

Run  and  calibrate 

Scale  Switch 

5,  10,  20  m/sec 

Aperture  Switch 

Small,  medium,  large 

Time  Constant  Switch 

1 , 10,  100  sec 

Field  of  View 

0,6  degree 

Optics 

Twin  apertures  with 

Fresnel  lenses 

Detectors 

United  Detector 

Technology  pin  spot  2D 
Photodiodes  (pair) 

Size 

25  X 30  X 45  cm 

Power  Requirements 

120  V AC 

Transmitter 

Light  Source 

Quartz-iodine  55  W 12  V DC 
at  focus  of  Fresnel  lens 

Power  Requirements 

120  V AC 
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Figure  2.  Saturation  Resistant  Crosswind  Sensor  Block  Diagram. 


Figure  3.  Saturation  Resistant  Crosswind  Sensor  Transmitter  (oblique  view). 


The  operation  of  the  SRCS  is  based  on  detection  of  the  scintillation 
patterns  produced  by  thermal  gradients  and  transported  by  the  wind. 

As  each  one  of  the  photodiodes  detects  the  incoming  signal,^  the  fluctua 
tions  in  irradiance  are  combined  by  the  analyzer  to  obtain  a covariance 
function.  The  slope  of  this  curve,  at  zero  delay  time,  is  proportional 
to  the  time  needed  for  a particular  scintillation  pattern  to  travel  from 
one  detector  to  the  other.  Since  the  distances  between  photodiodes  are 
constant,  the  analyzer  then  derives  the  wind  velocity  from  the  slope  of 
the  covariance  function.  Figure  2,  a block  diagram  of  the  SRCS,  illus- 
trates the  electronic  circuitry  required. 


TEST  SUPPORT 
Biggs  Optical  Range 

The  BOR  is  located  approximately  400  m northwest  of  the  main  runway  at 
Biggs  Army  Airfield,  Fort  Bliss,  Texas.  The  range  consists  of  an  instru- 
mented path  at  a heading  of  49  degrees  from  True  North  and  2064  m long. 
Two  3.5-m  towers  are  located  at  the  end  points  of  the  path  with  a 3-m 
tower  aligned  at  the  500-m  point  of  the  path.  These  towers  provide 
solid  test  beds  for  electro-optical  instrumentation.  A linear  array  of 
3-m  high  anemometers  parallel  to  the  path  is  offset  3 m to  the  southeast 
of  the  optical  path.  The  anemometers  are  oriented  to  measure  northwest- 
southeast  winds,  i.e.,  "cross"  winds  to  the  optical  path.  This  array 
consists  of  21  anemometers  spaced  25  m apart  for  the  first  500  m path 
length  and  15  anemometers  spaced  100  m apart  for  the  remaining  1500  m 
of  path  length.  Figure  4 shows  the  BOR  with  the  SRCS  in  the  foreground. 

The  instrumentation  path  is  specifically  designed  to  test  optical  wind 
measurement  systems  because  the  surrounding  terrain  features  are  flat 
and  the  optical  path  has  been  cleared  of  natural  vegetation  to  minimize 
wind  flow  field  characteristics. 

Meteorological  Optical  Measuring  System 

The  Meteorological  Optical  Measuring  System  (MOMS)  (figure  5)  is  a 
mobile,  self-contained  data  collection  and  reduction  system  containing 
analog  and  digital  subsystems  specifically  engineered  for  the  measurement 
and  recording  of  atmospheric  meteorological  data.  The  system  uses  an 
HP  2100  computer  system  as  a controller  and  is  managed  by  an  in-house 
developed  program  that  samples  the  various  sensors  at  present  rates. 


®G.  R.  Ochs,  S.  F.  Clifford,  and  Ting-i  Wang,  1976,  "Laser  Wind  sensing: 
the  effects  of  saturation  of  scintillation,"  Appl  Opt,  15,  No.  2,  403-408 
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Figure  5.  Meteorological  Optical  Measuring  System. 


stores  these  data,  and  then  reduces  and  analyzes  these  data  according  to 
the  program.  Output  format  capabilities  are  raw  scatter  graphs,  time 
averaged  plots,  printer,  limited  strip  chart,  and  digital  tape. 

During  these  tests,  analog  wind  data  from  the  anemometer  array  and  the 
SRCS  output  were  recorded  on  digital  tape.  Other  meteorological  data 
simultaneously  recorded  were  atmospheric  pressure,  temperature,  refrac- 
tive index  structure  coefficient,  and  dew  point. 

As  a part  of  the  data  collection  and  analysis  effort,  data  reduction  was 
conducted  both  on  and  off  line  by  the  FORTRAN  program  shown  in  appendix 
D.  The  primary  results  provided  were  scatter  plots  and  weighting  func- 
tion diagrams  of  the  SRCS  versus  the  anemometer  array. 


Remote  Sensing  Van 

The  Remote  Sensing  Van  (RSV)  is  a 5-ton,  6 by  6,  M820  expandable  van 
which  contains  inherent  prime-mover  mobility  and  provides  test-bed 
facilities.  The  van  "folds"  to  standard  van  width  for  transport  and 
expands  to  4.3  m for  in  situ  operation.  The  RSV  is  a stable  platform 
for  optical  equipment  tests  and  provides  test-bed  facilities  by  housing 
test  equipment  and  ancillary  dedicated  test  support  equipment.  An 
environmental  isolation  screen  with  two  30-  by  45-cm  integral  glass 
plates  has  been  fabricated  for  use  so  that  the  rear  doors  can  be  opened 
for  optics  1 ine-of-sight  test  capability  while  test  environmental  condi- 
tions are  retained  inside  the  RSV.  Figure  6 shows  a "downrange"  view 
of  the  RSV  in  operating  configuration. 


TEST  DESCRIPTION,  CONDUCT,  AND  PROCEDURES 
General 

The  evaluation  mission  was  twofold:  (1)  determining  the  accuracy  and 
weighting  function  of  the  SRCS  and  (2)  testing  SRCS  operational  charac- 
teristics for  effects  due  to  vibration  and  weather  conditions  (i.e., 
rain,  overcast). 

Range  operational  setup  of  the  SRCS  involved  placing  the  transmitter  at 
one  end  of  the  path  and  the  receiver  at  the  other  end.  To  insure  opera- 
tion, the  transmitter  and  receiver  should  be  focused  accurately  toward 
each  other.  The  transmitter  could  be  aligned  easily  and  rapidly  by 
placing  a 7.5-cm  retroref lector  on  the  receiver.  After  the  transmitter 
was  aligned,  the  receiver  could  be  successfully  aligned  by  observing 
electrical  signals  from  the  "test"  outputs  while  slowly  slewing  the 
receiver  in  azimuth  and  elevation  until  the  strongest  scintillation 
signals  were  observed  from  both  channels.  When  the  alignment  procedure 
was  completed,  both  receiver  and  transmitter  units  were  secured  to  avoid 
any  deterioration  of  the  signal -to-noise  ratio  (SNR)  and  any  "false" 
signal  generation  due  to  vibration. 
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Figure  6.  Remote  Sensing  Van. 


Ranges 


The  SRCS  was  tested  at  500-m  and  2000-m  path  lengths.  SNR  at  500  m was 
relatively  high  and  proper  lockon  was  readily  obtained.  Extending  the 
operating  range  to  2000  in  decreased  the  SNR,  but  by  proper  transmitter- 
receiver  alignment,  lockon  was  achieved. 


Weather  Conditions 

Because  the  SRCS  is  an  active  system,  its  operation  is  not  limited  to 
daytime  use.  In  fact,  the  ope'  ating  range  increases  at  night  because  of 
the  absence  of  ambient  light  (considered  extraneous)  that  tends  to  degrade 
the  SNR.  Operation  on  clear  or  overcast  days  was  reliable;  however,  the 
operation  became  erratic  during  rain  and  ceased  during  heavy  rain. 

A detailed  daily  summary  of  weather  conditions  existing  during  the  test 
period  is  shown  in  appendix  C.  These  data  are  from  the  National  Weather 
Service  located  at  the  El  Paso  International  Airport  approximately  6 km 
from  BOR.  A synoptic  weather  summary  for  the  surface  and  500  mb  altitude 
is  also  shown. 

Atmospheric  conditions  of  thermal  turbulence  ranged  from  very  low  or  non- 
existent (during  periods  of  rainfall)  to  values  large  enough  to  cause 
saturation  of  other  crosswind  systems. 


DATA  COLLECTION  AND  RESULTS 
Mathematical  Background 

In  statistical  bivariate  analysis,  a scatter  plot  is  useful  for  the  evalu- 
ation of  measured  experimental  data.  During  these  tests,  scatter  plots 
were  generated  because  the  resultant  comparison  values  were  important  for 
determination  of  accuracy. 

For  the  scatter  plots  employed,  the  average  path  crosswinds  measured  by 
groupings  of  weighted  anemometer  outputs  from  the  array  are  plotted  on 
the  abscissa,  and  the  measured  average  output  from  the  SRCS  is  plotted 
on  the  ordinate.  These  two  values  were  plotted  as  an  ordered  pair.  To 
make  this  comparison  applicable,  a large  number  of  these  single  sets  of 
values  had  to  be  compared.  This  comparison  should  result  in  a statis- 
tically sound  conclusion.  The  FORTRAN  IV  program  (appendix  D)  was 
developed  to  sample  and  plot  the  experimental  data. 

The  usefulness  and  simplicity  of  the  scatter  diagram  is  readily  shown  in 
figure  7.  The  first  scatter  plot  (plot  A)  shows  a straight  line  with  a 
45-degree  slope  and  passing  through  the  origin.  In  this  case,  a one-to- 
one  correspondence  (complete  agreement)  exists  between  the  contributing 
systems.  Plot  B differs  in  that  the  slope  of  the  line  is  no  longer  45 
degrees.  This  case  indicates  that  the  ordinate  values  have  to  be 
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AWA  VALUE 
=SRCS  VALUE 


Typical  scatter  plots.  Windspeed  values  measured  by  the  SRCS  are  plotted  versus 
values  measured  by  a reference,  in  this  case,  the  analog  wind  averager. 


adjusted  by  including  a constant  multiplicative  factor  (m).  Plot  C shows 
the  resultant  line  no  longer  passing  through  the  origin.  This  signifies 
the  ordinate  values  have  to  be  further  adjusted  by  including  an  offset 
value  Yq. 

Since  experimental  data  rarely  yield  complete  correspondence  between 
test  system  and  "base"  systems,  data  obtained  were  linearized  by  employ- 
ing the  least  squares  fit  method.  The  result  was  then  a line  represented 
as 

Y = A + A,X, 
o 1 

where^ 

A = - (SXXZXY) 

° NCix^)  - (ZX)^ 


and 


NgXY)  - (ZXXZY) 
NCZX^)  - (EX)^ 


N = number  of  samples. 

The  method  of  least  squares  allows  extreme  values  to  weigh  too  heavily  on 
the  result;  therefore,  when  this  technique  was  used  to  evaluate  the  SRCS, 
as  well  as  other  test  systems,  care  was  taken  to  investigate  extremely 
large  differences  (on  actual  measurements)  to  insure  that  they  were  legit- 
imate and  not  caused  by  mechanical  or  electrical  malfunctions  of  one  of 
the  systems.  Subjective  caution  also  allowed  for  the  fact  that  the 
mechanical  anemometers  have  a certain  amount  of  inertia  which  normally 
results  in  an  erroneous  output  when  the  wind  velocity  is  below  the 
threshold  value  of  the  anemometers.  Therefore,  it  is  necessary  to  con- 
sider the  values  close  to  zero  carefully  and  prudently.  The  anemometers 
used  minimized  this  anomaly,  however,  since  they  are  research  quality 
propeller  anemometers  with  a threshold  value  of  0.2  - 0.3  m/sec. 


^Anthanasios  Papoulis,  1965,  Probability,  Random  Variables,  and  Stochastic 
Processes,  McGraw-Hill  Book  Company 
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Scatter  plots  can  be  generated  by  using  either  straight  average  or 
weighted  average  values  from  the  analog  wind  averager  (AWA)  as  the 
abscissa  point  input.  First,  plots  were  generated  by  using  straight 
average  values;  in  later  tests,  after  the  weighting  function  of  the 
SRCS  had  been  determined,  scatter  plots  with  weighted  average  values 
from  the  AWA  were  generated. 

Both  types  of  plots  are  illustrated  in  appendix  A;  however,  more 
emphasis  is  given  to  the  weighted  value  plots. 

Before  these  weighted  values  could  be  used,  the  SRCS  weighting  function 
had  to  be  determined.  The  weighting  functions  were  computed  by  con- 
sidering different  groups  of  anemometer  outputs  as  a least  squares  basis  set; 
that  is,  the  SRCS  wind  measurements,  denoted  by  W^,  are  represented  as 

a linear  combination  of  m different  groups  of  anemometer  outputs,  W^. , 
so  that 


m 

w = . a.W.  , 

S 1=1  1 l’ 


where  a^.  are  the  correlation  coefficients  determined  by  the  m^*^  order 

least  squares  analysis.  Various  sets  of  coefficients,  a.,  were  obtained 

by  employing  different  groupings  of  anemometers;  and  since  the  range  is 
2 km  long,  it  was  possible  to  obtain  a weighting  function  out  to  that 
distance. 


Results  of  Data  Analysis 

The  weighting  function  and  scatter  plots  obtained  from  the  collected 
data  are  compiled  in  appendix  A for  easy  reference.  Figure  8 is  a repre- 
sentative plot  that  exemplifies  the  method  of  data  presentation.  For 
this  case,  the  SRCS  data  closely  followed  the  AWA  data,  as  can  be  seen 
by  the  dashed  line  which  approaches  quite  well  the  solid  (45-degree)  line. 
The  weighting  function  graph  illustrates  the  preponderance  of  weighting 
from  250  m to  400  m with  a maximum  at  approximately  350  m.  For  the  sake 
of  completeness,  figure  8 is  replicated  and  included  with  the  rest  of  the 
reduced  data  in  appendix  A.  In  summary,  for  a path  length  of  500  m the 
values  of  the  SRCS  are  related  to  those  of  the  AWA  if  the  following 
formula  Y = 0.85X  - 10.02  and  the  weighting  function  show  that  cross- 
winds  located  from  250  to  400  m are  weighed  most  heavily.  For  the 
2000-m  path,  the  previous  results  were  shifted  somewhat  to  yield 
Y = 1.57X  - 0.43  and  the  path  length  weighed  most  heavily  is  from 
1000  to  1700  m. 

Appendix  B contains  the  comparison  graphs  of  the  SRCS  versus  the  AWA  as 
a function  of  time.  Comparison  of  these  values  yielded  a measure  of 
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accuracy.  Calculations  show  that  the  crosswind  values  measured  by  the 
SRCS  were  within  8 percent  of  those  measured  by  the  AWA  at  least  81  percent 
of  the  time. 

The  program  included  in  appendix  D was  used  to  collect,  reduce,  and  plot 
the  data  presented  in  this  report  and  is  not  discussed  per  se  in  this 
report. 

Although  the  weighting  function  derived  for  the  SRCS  with  a 500-m  range 
is  not  exactly  the  same  as  the  one  derived  with  a 2000-m  range,  values 
calculated  from  them  are  within  15  percent  of  each  other.  To  determine 
the  exact  cause  of  the  discrepancy,  more  data  than  what  were  collected 
during  the  test  period  would  have  to  be  gathered  and  analyzed  statistically 

CONCLUSIONS 

Under  predominant  area  weather  conditions  (no  rain,  hail,  or  snow),  the 
SRCS  measured  crosswind  averages  within  8 percent  of  those  measured  by 
the  AWA.  However,  under  moderate  to  severe  adverse  weather  conditions, 
the  SRCS  was  not  accurate  or  reliable  in  operation,  due  to  a decrease  of 
SNR  below  its  operating  threshold.  Both  units  of  the  SRCS  should  be 
solidly  attached  to  a firm  foundation  to  eliminate  any  vibration,  since 
such  vibration  will  be  interpreted  as  crosswind  signal  information  and 
thus  yield  erroneous  results. 

The  calculated  SRCS  weighting  functions  indicate  that  crosswinds  closer 
than  1/5  of  the  path  from  the  receiver  or  closer  than  VIO  of  the  path 
from  the  transmitter  affect  its  output  reading  minimally,  while  cross- 
winds  at  the  central  section  of  the  path  affect  the  output  predominately. 

The  SRCS  transmitter  can  easily  be  aligned  by  using  a retroref lector  at 
the  receiver  and  a boresighted  telescope  at  the  transmitter  end.  The 
receiver  can  then  be  aligned  by  observing  the  signals  out  of  the  "test" 
connectors  and  minutely  slewing  the  receiver  slowly  until  maximum  scin- 
tillation readings  are  obtained  at  both  connectors. 

The  SRCS  can  aid  in  the  characterization  of  the  atmosphere,  serve  as  a 
reference  to  evaluate  future  wind  measurement  systems,  and  help  advance 
research  on  wind  measurement  techniques. 

The  use  of  the  SRCS  is  particularly  applicable  for  crosswind  measurements 
during  conditions  of  high  thermal  turbulence  or  over  long  path  lengths 
where  other  crosswind  measurement  systems  have  exhibited  optical  saturation 

Two  weighting  functions,  instead  of  one,  were  derived  for  the  SRCS  to  have 
a higher  accuracy  of  values  read.  However,  if  only  one  is  desired,  an 
average  weighting  function  could  be  obtained,  with  a corresponding  decrease 
in  accuracy. 
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APPENDIX  B.  SRCS  WIND  MEASUREMENT  COMPARISON  PLOTS 
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Figure  B-3.  SRCS  wind  comparison  plot,  13  Feb  78,  0.5  km  range,  1115-1215  hours 


TIME 

Figure  B-4.  SRCS  wind  comparison  plot,  13  Feb  78,  0.5  km  i-nge,  1215-1315  hours. 
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Figure  B-5.  SRCS  wind  comparison  plot,  13  Feb  78,  0.5  km  range,  1315-1325  hours. 
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Figure  B-7.  SRCS  wind  comparison  plot,  28  Feb  78,  0.5  km  range.  1335-1435  hours. 
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Figure  B-8.  SRCS  wind  comparison  plot,  28  Feb  78,  0.5  km  range,  1435-1535  hours. 


TIME 

Figure  B-9.  SRCS  wind  comparison  plot,  28  Feb  78,  0.5  km  range,  1535-1555  hours. 


APPENDIX  D.  FORTRAN  IV  DATA  PLOT  PROGRAMS 


JULTHC  T=>.'0v04  is  oh  CPOt'vvt:  USING  0O020  BLKS  fi;=014i& 

1 FTN4.L 

o-:K'2  PPOGRhM  WLTHC,3 

<>-;:'f)4  c WLTMC  IS  USED  TO  LOmD  COMMON  WITH  THE  DESIRED 

D005  C FARHMETERS  FOR  MAG  TAPE  RNHLYSIS  WITH  WLTMG 


0 0 0 7 
0 0(>8 
000  8 
0 0 1 0 
0 011 
001  2 
0013 
0014 
0015 

0 0 1 s 

0017 
0 0 1 8 
0 0 1 9 
002  0 
0 021 
0 022 
0023 
002  4 
002  5 
0026 


COMMON  F , S( 22 > . I T( 6 ) 

COMMON  A7«:  4 , 4 ) , B7(  4 XP7i:  4 NH7  , NS7(  4 ) , SB7.  BA7(  4 > 
COMMON  A6(5..  5 ),B6<  5 XP6<5  ),  NH6..NS6(  5).  SB6..  886(5  I 
COMMON  85(5.5:>,B5<5).XP5(5>..N85.NS5<5>,SB5..  B85(5) 
COMMON  84(6.6  >,B4(  6 ).  XP4(6  N84..NS4<  6 > .■  SB4.  B84(6  > 

COMMON  8 3(8. 8 ),e3( 8 ). XP3(e  ). N83. NS3( 8 >. SB3. B8  3<  8 ) 
COMMON  82(  1 1 . 1 1 > B2(  1 1 >.  XP2(  1 1 ) .•  NA2.  NS2(  1 1 > . SB2.  e82(  1 
COMMON  P, NB.MN, IS. MS. I8RRY 
COMMON  IWl.  U12, IW3.  IW4.  I W5 . IWZ 
COMMON  WT(21> 

COMMON  XC. YC. SN.SX. SY.SXX.SYX.DA. DS 
DIMENSION  IP(5> 

CALL  RHF'AR(IP> 

LUl  = IP( 1 ) 

LU2=  IP( 2 > 

LU3  = IP(3  > 

LP=IP(  4 1' 

IFv  LUl  . EQ  0 )LU1  = 1 
IF( LU2  ED  O )LU2=  1 
IF( LU3 . EQ  .0  )LU3=  1 


0027 
0 02  8 C 
0 02  9 
0 0 3 0 
0 0 3 1 
0 03  2 
003  3 
0 0 3 4 
003  5 
0 03  6 
0 03  7 
0038 
0 03  9 
004  0 


IF( LP  EQ . 0>LP=1 
2ER0  OUT  ARRAYS 
DO  1 1=1.11 
DO  1 0=1.11 
A2(  I , 0 > = 0 . 0 
B2(  I > = 0 0 
XP2<  n=o . 0 
NS2(  I ) = 0 
EA2(I>=0  0 

I F(  I GT  8 OR  0 . GT  8 )G0  TO  1 

A3(  I , 0 ) = 0 0 

B3(  I ) = 0 0 

XP3(  I > = 0 . 0 

NS3(  I •>  = 0 


004  1 
0042 
0 04  3 
0 04  4 
004  5 
0 04  6 
0 04  7 
0048 
004  9 
0050 
0051 
005  2 
0053 
005  4 
0055 
005  6 
005  7 
0058 


BA3(  I > = 0 . 0 

IF(  I GT . 6 OR  J . GT  6 )G0  TO  1 
A4<  I . J ) = 0 . 0 
B4<  I ) = 0 . 0 
XP4C  n = o 0 


MS4(  I ) = 0 


BA4<  I ) = 0 . 0 


IF(I.GT.50R.J.GT.5)G0  TO  1 

85(  I , 0 > = 0 . 0 

B5(  I ) = 0 . 0 

XP5C  I > = 0 . 0 

NS5(  I ) = 0 


BA5(  I >=0 . 0 
86(  I . J ) = 0 . 0 
B6<  I ) = 0 . 0 
XP6(  I )=0  . 0 
HS6(  1 ) = 0 
BA6(  I ) = 0 0 


mis  PiGiL  lb  4„ALlTy 

PROW  OOr'Y  r 10  ODC 
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•')  05  9 

I F«  I GT  . 4 OR  J . GT  4 )G0  TO  1 

0 0 fe  0 

A7(  ! . J ) = ''  0 

•v'Vb  1 

e7<'  I > = 0 0 

0 06  2 

X P 7 (.  ! > = 0 . 0 

0 06  3 

MS7(;  I ) = 0 

<>  0 6 4 

8 371.  10  = 0.0 

A ^ 5 

j 

CONTINUE 

0 0 b b 

NA7  = 0 

0 0 6 7 

887=0  0 

0 0 6 S 

Nh6  = 0 

0 06  9 

886=0  0 

0 0 7 0 

NA5  = 0 

0 0 7 1 

885  = 0 . 0 

()  Y>  7 2 

NA4  = 0 

0 0 7 3 

884  = 0 . 0 

0 0 ? 4 

NA3  = 0 

0 0 7 5 

883=0 . 0 

0 0 7 6 

HA2  = 0 

0 0 7 7 

8 82  = 0.0 

0 07S 

c 

ZERO  SUMS  FOR  ULTMl 

007  9 

P = 0 . 0 

0 0 e 0 

DA  = 0 . 0 

008  ! 

DS  = 0 . 0 

008  2 

SN  = 0 . 0 

0 08  3 

SX  = 0 . 0 

0 08  4 

SY  = 0 . 0 

008  5 

SYX  = 0 . 0 

0 08  6 

SXX  = 0 . 0 

0 08  7 

c 

GET  CALIBRATION  AND  WEIGHTS  FOR  ULTMl. 

0 08  8 

IF( LP. NE  1 )WRITE(LP, 199 ) 

0 08  9 

199 

FORMAT< IHl. 2X.  “CALIBRATION  FACTORS  FOR  REAL 

TIME  ANALYSIS" 

0 09  0 

+1H0, 2X, "WEIGHTS  FOR  ANEMOMETERS  ARE"/) 

0 09  1 

IF(  LUl  .EQ  . 1 )URITE(  1,99) 

0092 

99 

FORMAT< " INPUT  WEIGHTS"/) 

009  3 

DO  2 1=1.21 

0 09  4 

I F<  LUl  . EQ . 1 )WR  ITE(  1 , 98  ) I 

0095 

READ<  LUl . * )WT<  I ) 

0 09  6 

I Ft  LP. NE . 1 )URITE<LP, 198  TI- UT(  I ) 

0 09  7 

2 

CONTINUE 

0098 

98 

FOR.MAT<  I 2 , 3X.  " ) 

0 09  9 

193 

FORMAT< 1 H . 12. FI  0 . 5 ) 

0 10  0 

IF( LU2  EQ . 1 )URITE(  1.97) 

0 10  1 

97 

FORMAT< " INPUT  XCAL.  Y C A L " . 3K . " _ “ ) 

010  2 

READ<LU2. *)XC. YC 

0103 

IF<  LP . NE  . 1 )WRI TE<  LP , 1 97  )XC , YC 

0104 

197 

FORMAT< IHO. 2X. "XCAL. YCAL  ARE " , F 1 0 . 7 . “ . " . F 1 0 

. 7 ) 

0 105 

IF(  LU2  EQ . 1 )URITE(  1.96) 

0 106 

96 

FORMATC " INPUT  « OF  SE NS OR " , 3 X . " _ " ) 

0 lo7 

READ<LU2.  ■♦)IS 

0108 

IFCLP.NE  1) WRITE <LP,  196) IS 

0 109 

196 

FORMAT< IHO. 2X. "SENSOR  #",i4) 

0 110 

IF<  LU2  EQ . 1 )URITE(  1.95) 

0111 

95 

FORMAT( " INPUT  TIME  INTERVAL  IN  MINUTES  FOR 

ULTM3" .3X. ) 

0112 

READ(LU2.*)MN 

0113 

IF( LP. NE  1 )URITE<LP. 195 )NN 

0 114 

195 

FORMAT< IHO. 2X.  "LEAST  SQUARES  READOUT  EVERY" 

.13,"  MINUTES") 

0115 

C 

DETERMINE  IF  . SK  OR  2K  WANTED. 

0116 

IF( LU2 .EQ . 1 )WRITE< 1.93) 

0117 

93 

FORMAT( "ENTER  0 FOR  . 5K  OR  2 FOR  2K",3X."." 

) 

0118 

READ<LU2. * )IARRY 

mis  PAGE  IS  BiSI  QUALITY  PRAftlWiABIil 
i IfJfliHsHJC)  ro  DDC 


0 119 

0120 

0121 

193 

0 122 

200 

0123 

0 124 

92 

0 125 

0126 

0127 

192 

0128 

91 

0129 

90 

0130 

89 

0 131 

88 

0132 

87 

0133 

191 

0 134 

190 

0 135 

189 

0 136 

188 

0 137 

187 

0138 

0139 

0 140 

0141 

0142 

0 143 

0 144 

0 145 

0 146 

0 147 

0 148 

0 149 

0150 

0 15  1 

0 152 

0 153 

0 154 

35 

0155 

0 156 

0 157 

185 

0158 

0 159 

34 

0 160 

0161 

0 162 

184 

0 163 

0164 

83 

0165 

0166 

0 167 

IF(IARRY.EQ  0)WRITE<LP.193) 

IF(  I ARRY  . EQ  2 )lilRITE<LP.  200  ) ! 

F0RMAT< IHO, 2X, "ANALYSIS  FOR  THE  1/2  K ARRAY") 

FORtlATC  IHO.  2X.  "ANALYSIS  FOR  THE  2 K ARRAY")  i 

IF( LU3  EG . 1 )«RITE(  1,92) 

FORMATC " INPUT  0 OR  1 FOR  SUBS  DESIRED"/ 

+"0R  TO  DESIGNATE  OTHER  PARAMETERS") 

IF(  LP . NE . 1 )WR I TE<  LP  < 192  ) 

FORMAT< 1 HO, 2X , "SUBS  DESIRED  OR  OTHER  PARAMETERS") 

FORMAT< " lUl ", 3X, ) ; 

FORMAT< " IW2", 3X,  "_"  ) 

FORMAT< " IW3", 3X,  “_"  ) ! 

FORMATC " IW4", 3X,  "_"  ) 

FORMAT< " IU5", 3X,  "_"  ) 

FORMATUH  , 2X  , " I « 1 " , I 5 ) 

FORMATC 1 H , 2X, "IW2" , 15) 

FORMATCIH  , 2X , " I W3 " , I 5 ) 

FORMAT<  IH  , 2X,  “HJ4"  , 15) 

F0RMAT<1H  , 2X, "IU5" , 15) 

IF<  LU3 . EG  . 1 )URITE<  1 , 91  ) 

READ<LU3, *>IW1 

IF( LP . NE . 1 )WRITE<LP, 191  )IW1  ; 

IF( LU3  EG . 1 )WRITE(  1 , 90)  ’ 

READ<LU3,  f )IU2 

IF( LP . NE . 1 )«R I TE< LP , 1 90  ) IW2  I 

IF( LU3  EG . 1 )WRITE(  1,89) 

REA0<LU3,*)IW3 

IF( LP . NE . 1 )«R I TE< LP , 1 89  ) IU3  . 

IF( LU3 . EG . 1 )MRITE<  1,88) 

READCLUS, *)IW4 

IF(  LP . NE , 1 )WRI TEC LP , 1 88 ) IW4 

IF(LU3  EG. 1 )WRITE(  1 ,87)  i 

REA0CLU3,H.)iy5  ' 

IFC LP . NE , 1 )WRI TEC LP, 1 87  )IU5 
WRI  TEC  1,85) 

FORMATC " ENTER  FILE  # ON  TAPE " , 3X,  " _"  ) i 

READ<1,*)IWZ  I 

IFC  LP , NE . 1 )UR I TEC  LP , 1 85 ) IWZ 
FORMATC IHO, 2X, "TAPE  FILE  #",I2) 

WRI  TEC  1,84) 

FORMATC  " I NPUT  TIME  INTERVAL  IN  SE C ONDS " , 3 X , " _ " ) 

READC 1 , * )MS 

IFC  LP . NE  1 )WR ITEC  LP , 184  )MS 

FORMATC IHO, 2X, "AVERAGING  TIME  IS", 13,"  SECONDS") 

WRI TEC  1,83) 

FORMATC/"THAT'S  ALL,  THANKS") 

STOP 

END 

END* 


XBIS  PAM  i-  ' *.  • 
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tULTHG  T-00004  IS  OH  CR00002  USING  00014  BLKS  R>0114 


0001 
0002 
0003 
0 00  4 
0005 
0006 
000  7 
0008 
0009 
0010 
001  1 
0012 
0013 
0014 
0015 
0016 
001  7 
0018 
001  9 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0 03  0 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 
0048 
0049 
0050 
0 05  1 
0052 
005  3 
0 05  4 
005  5 
0056 
0057 
005  8 


FTN4,L 


PR0GR8H  yLTHG,3 


COHNON  F/ S( 22  >. 1 T( 6 ) 

C OH  HON  R7<4,4  ),B7(  4 ),  XP7<4  ),  NA7  , HS7<  4 ) , SB7.  BA7<4  ) 

C OH HON  86(5.5  ).66<  5 >. XP6( 5 > . NA6 . NS6( 5 ) . SB6. B86< 5 ) 

COHHON  85(5 . 5 >.B5( 5 >. XP5(5  ). N85. NS5( 5 > . S65. 685( 5 > 

COHHOH  84(6.6)>B4(6).XP4(6  ) . N84 . NS 4( 6 > . SB 4 . 684(6  > 

COHHON  83( 8 . 8 ) >B3( 8 >. XP3(8  >. N83. HS3( e>. SB3. B83( 8 ) 

COHHON  82( 1 1. 1 1 >.B2( 1 1 ). XP2( 1 1 >.NA2. NS2( 1 1 >. SB2. B82( 1 1 > 
COHHON  P. NB.HN. IS. HS. I8RRY 
COHHON  IUl.iy2.IH3.IU4.IU5.IW2 
COHNON  WT(21  > 

COHHON  XC.YC.SN.SX.SY.SXX.SYX.DA.DS 
D I HENS  ION  ITN(5>.IDT(94>.ID8T8(  100> 

DIHENSION  ST(37) 

DIHENSION  NU4( 3).H0N( 2> 

EQUIVALENCE  ( I TH( 1 ) . I D8T8(  1 ) ) . ( I Y . ID8T 8( 6 ) ) . 

♦ ( IOT(  1 ).  ID8T8(  7)) 

DATA  NU4/2HUL . 2HTH . 2H4  / 

CONTINUE 

CALL  EXEC(3.611B) 

CALL  EXEC( 13.9. ISTAT) 

1STT  = I AHD( ISTAT. IB  ) 

IF( ISTT  NE  0)G0  TO  1 
ISTT=IAHD( ISTAT. 200B) 

IF( ISTT . EQ  0)G0  TO  2 
CALL  EXEC(3.311B) 

CONTINUE 

CALL  EXEC( 1 . 1 1 IB. IDATA. 100 ) 

CALL  EXEC(3.211B) 

SEC=HS 

IF( HS . EQ . 0)SEC=0  5 
IDAY=ITH( 5) 

LIT  = ITM( 3 ) 

CALL  DATE( IDAY.HON.  lY  ) 

URITE(6.99)IS.SEC.  IU3.ITH(4).  ITH(3).  ITH(2).  IDAY.HON.  I Y 


WITH" .F5  . 1 


99  FORHAT( IHl. "ANALYSIS  OF  CH  •".13.". 

*1H  ."WITH  SLIDE  FACTOR  OF". 13/ 

♦ IN  ."FOR  DATA  BEG I NNI NG " . I 3 . " : " . I 2 . 

CALL  EXEC( 1 1.  ITH. lY  ) 

IOAY=ITH( 5) 

CALL  DATEdDAY.HON.  lY  ) 

URI TE( 6. 98)  ITH( 4 ).  ITH( 3 ) . I TH( 2>. IDAY.HON. lY 
98  FORHATdH  ."ANALYSIS  ST  ARTED  " . I 3.  * ! " , I 2 . " s * . 12 
IF(  lUl  EQ  . 0 )G0  TO  3 
CALL  ULTHS 

3 CONTINUE 

DO  67  IRPT=1. 32767 
F=0 . 0 

DO  4 1=1.37 

4 ST(I)  = 0.0 

DO  6 1=1 . 200 

CALL  EXEC( 1 . 1 1 IB. IDATA. 100 ) 

CALL  EXEC(13.9.ISTAT) 

ISTAT=IAND(  ISTAT. 200B  ) 

IF( ISTAT  NE .0)G0  TO  68 
F=F+1 


SEC  AVG."/ 


.12."  OH". 13. IK. 2A2. 14) 


ON*.  13.  IK. 2A2.  I 4 
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0 C)  5 9 

0 0 6 0 

5 

0 06  1 

0 06  2 

0 06  3 

C>06  4 

0 0 6 5 

0 0 6 6 

0 06  7 

6 

0 C>  6 8 

7 

0C>69 

0070 

8 

0071 

0072 

0 0 7 3 

007  4 

007  5 

9 

007  6 

C 

0 07  7 

0078 

0 07  9 

008  0 

0 081 

10 

008  2 

008  3 

008  4 

1 1 

0 08  5 

0 08  6 

12 

0 08  7 

OoSS 

13 

008  9 

14 

009  0 

009  1 

0 

0 09  2 

009  3 

0 094 

0 0 9 5 

15 

C>  *0  9 6 

0097 

009  3 

0099 

0 1 0 0 

0101 

0 102 

0 ! 0 3 

67 

0 1 0 4 

f P 

0 105 

0106 

0 107 

0 108 

0 1 ’.  9 

. 16 

0 110 

1 7 

0 111 

0112 

0 113 

18 

0 ! M 

0 115 

0116 

6 9 

0117 

0 118 

DO  5 J=1 , 36 

ST( 3 )=ST<  0 ) + FLOhT(  IDT( J > ) 
ST(3?)  = ST<37)  + FLOflT(IDT(  IS)> 
IFC  MS. EQ . 0)G0  TO  7 
ISTftT  = MOD< I TM<  2 ) ,nS  ) 

IF(  ISTFiT  NE  .0  >G0  TO  6 

IFC  I TN< 1 ) . GE . 50  )G0  TO  6 

GO  TO  7 

CONTINUE 

CONTINUE 

00  8 J=1 , 6 

I T(  J )=IDflTf(f  J ) 

IFC  IT( 3)  . EO  .LIT)GO  TO  9 
ISWTM=ITC  3>  + ITC4  )*1000e 
CALL  PSSUCISWTM) 

L IT=^  ITC  3 ) 

CONTINUE 

CHECK  FOR  . 5K  OR  2K  RANGE 
IFC  lARRV  LE . 1 )G0  TO  12 
K = 0 

DO  10  1 = 1 ,21..  4 


K=K+  1 

SCK  ) = STC  I )*XC/F 
DO  11  1=22,36 
K = K+  1 

SCK  ) = STC  I )*XC7F 
GO  TO  14 
CONTINUE 
DO  13  1=1,21 
SCI  ) = STC  I >*^;CC,''F 
CONT INUE 

SC  22  >=STC  37  >*YC/F 
SUM  ALL  POINTS 
P = P + F 

IFC  IWl  EO .0  )G0  TO  15 

CALL  WLTMl 

CONTINUE 

CALL  WLTM2 

TFCMN  ES.OJGO  TO  67 

MO  = MODC I TC  3 >, MN) 

IFC  MO . NE  0 )G0  TO  67 
IFC  ITC  2 > . NE  . 0 >GCI  TO  67 
I FC  I TC  1 ) GE  50  IGO  TO  67 
CALL  WLTMj 
CONI INUfc 
CONTINUE 

I K.  1 02  . L fc  . 1 )G0  TO  16 
LhLL  E XE CC 3, 2 1 IB  ) 

CAIL  EXE  CO,  141  IB) 

GO  TO  17 

CAI.  L EXECC  3 ,4  1 IE: 

CONT INUE 

IFC I H4 . HE . 1 )C0  TO  18 

CAIL  UlTHO 

CONT  INUE 

CALL  kllTM3 

CALL  EXECC 10, N«4  ) 

CONt INUE 

STOP 

END 


Oil?  END* 


SWLTMl  7=00004  IS  ON  CR00002  USING  OOOOt  BLKS  R=0043 


0 0 C-  1 
0 0 0 2 
0 0 0 i; 
0 0 0 4 

ooos 

0 0 V ii 
000  7 
0 0 0 8 
0 0 0 9 
0 0 1 0 
Ovi  1 
0 0 1 2 
0 0 1 3 
0 0 1 4 
0015 

V 0 1 t. 

0 01  7 
001  8 
0 0 1 9 
0 02  0 
0021 
002  2 
0023 
0024 
0 02  5 
0 026 
002  7 
0028 
0029 
0030 
0 03  1 
0 0 3 2 
0033 
0034 
003  5 
003  6 
0 05  7 
003  8 
0 03  9 
0040 
004  1 
004  2 
004  3 
0 04  4 
0045 


FI H4  , L 

SUBROUTINE  WLTMl 

0 » + :t  *■:»  ' T S ' I-  t * + * *■  (I  * + ++****■*♦♦  ♦*♦*•*  t *+*  + ♦*♦****♦♦♦♦♦♦*  If  ♦ If  ♦ * » 

C (fiLlt'll  IS  USED  TO  MftKE  THE  CORRELAIION  PLOT  OF  A SENSOR 

C VERSUS  A WEIGHTED  AVERAGE  OF  21  ANEMOMETERS  FROM  DATA 

C FROM  MAG  TAPE. 

Gum  MON  F,S<22>,n(6) 

COMMON  A7-:4  ..  4 B7<  4 ) / XP?<:4  ).  NA7  , NS7<  4 ),  SB?  / BA7<  4 ) 

COMMON  A6<5,5  >.B6<  5 ),  XP6<;5  ).•  NA6.NS6<  5>.  SB6..  BA6<5  ) 

COMMON  A5<  5.5) . B5(  5 >.  XPST  5 >.  NA5.NS5(  5).  SB5.  BA5<  5 ) 

COMMON  A4>;6  6 > .B4<  6 ).  XP4-'  6 >.•  NA4  . NS4(  6 > , SB4  , B A4  < 6 ) 

COMMON  ASC  8 . 8 ) .•  B3<  8 ).  XP3<  8 ),.  NA3..  NS3<  8 ).  SB3.  BA3-:  8 ) 

COM  MOM  A2<  1 1 ..  1 i ) . B2<  11  ).  X.P2(  1 1 ) ..  NA2.  NS2<  1 1 ) . SB2.  BA2<  1 1 ) 
COMMON  P . NEi ..  MN  ..  1 S..  MS..  I ARRY 
COMMON  IWl..  UI2.  IU3..  IW4.  IW5.  IW2 
COMMON  WT<21  ) 

COMMON  XC  ..  YC..  SN.  SX  . SY  . SXX.  SYX  . PA.  OS 
CAL  L PLTLU<  10  ) 

CALL  SFACTt;  15  . . 10  . ) 

Y = S<:  22  > 

X=-0 . 0 

DO  1 1=1.21 
1 X = X + S(  I )*WT< I ) 

SN=SN+ 1 . 0 
D = ( Y-X  ),'  1 0 . 

DA=DA+D 
DS=DS  + D*  D 
SX=SX+X 
SY=SY+Y 
SYX=SYX+Y+X 
SXX  = SXX+-X*X 
2 = ABS(  X ) 

IF( 2 . GE . 4 . 75)X=4 . 75*X/2 
X=X+5 . 0 
2 = ABS':  Y > 

IF<  2 . GE . 4 . 75)Y  = 4 . 75*Y/2 
Y=Y+5 . 0 

CALL  PL0T<X.Y.3) 

CALL  PL0T<X.Y.2) 

CALL  PL0T<X.y.3) 

69  CONTINUE 
RETURN 
END 
ENDt 
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SULTHS  T=00ft04  IS  ON  CR00002  USING  00021  BLKS  R=0205 


0001 
0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
0010 
001  1 
001  2 
0013 
0014 
0015 
001  6 
001  7 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
002  6 
0027 
0028 
0029 
0 03  0 
0031 
0032 
0033 
0034 
0035 
0 03  6 
0037 
0038 
0039 
004  0 
004  1 
004  2 
0043 
0044 
0 04  5 
0046 
004  7 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
005  7 
0058 


FI  N4  I 

SUBROUTINE  ULTN2 

C ULTN2  IS  USED  TO  LORD  THE  8RR8YS  FOR  THE  LEAST  SOUARES 

C CORRELATION  FIT  TO  NEIGHTING  FACTORS  FOR  DATA  FROfI  MAG 

C TARE 

COMMON  F . S<  22  >. I Tv  6 > 

COMMON  A7(4,4>,e7(4)<XP7<4),NA7,NS7<4).SB7,BA7<4> 

COMMON  A6<  5 ,5  ) .66<  5 ).  XP6<5  >,  MAS  ,MS6<  5>  < S66,  BAS';  5 T 
COMMON  A5t  5 <5 ),B5( 5 ), XP5<5 Nh5,NS5<  5 ) , SB 5 - B A5<  5 j 
COMMON  A4<6.6),B4(6;,XP4C6),NA4,NS4(6>,S6  4,BA4<6) 

COMMON  A3( 8-6  ),B3<  e T, XP3<  6 NA3- NS3( 8 ) , SB  3 - B A3 < 8 ) 

COMMON  A2(  1 1 - 1 n ,B2<  1 n - XP2<  1 1 ) ,NA2,  NS2<  1 n - SB2  - BA2<  1 1 ) 
COMMON  P,NB-MN,IS-MS- lARRY 
COMMON  lUl- IU2-1W3- IU4-  lyS-IUZ 
COMMON  WT<21) 

COMMON  XC-YC-SH-SX-SY,SXX-SVX-DA-DS 
DIMENSION  D( 1 1 )- NS<  1 1 >- XP(  1 1 > 

1F(  IU2  GT . 1 . )G0  TO  68 

HB  = 0 

NA7»0 

NA6  = 0 

NA5  = 0 

HA4  = 0 

NA3  = 0 

NA2  = 0 

DO  67  INOU^l- lUS 

GO  TO  < 1 - 2- 3- 4 , 5-6  ), INOU 

1 JCHT=4 
(<CNT  = 7 
GO  TO  7 

2 JCHT=5 
KCNT=6 
GO  TO  7 

3 JCHT=5 
KCNT=5 
GO  TO  7 

4 JCNT-6 
KCHT=4 
GO  TO  7 

5 JCNT=8 
KCNT=3 
GO  TO  7 

6 JCNT.=  11 
KCHT=2 

7 CONTINUE 
NAc  1 

DO  9 J=1 - JCNT 
NS<  HA)  = 0 
XP<  NAT  = 0 . 0 
DO  8 X=1 - KCNT 
L«( J-1 >*KCHT+K+I«3 


8 


I F( L . L T 1 OR  L . GT . 21  )C0  TO  6 
1FLAG=L 

NS( HA  ) = HS( NA  )+  1 

XP^  NA)  = XF<NA>  + FLOAT(L  ) 

CONTINUE 
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0 05  9 

1 FC  H3<:  Nft  > EQ  . 0 )G0  TO  9 

(>  0 t; 

MDEC  = 1!8<.  ho  ■' 

OOfc  1 

MFC  NO  ,'-XFC  NO  )/ODEC 

0 0 6 i 

NO=N0+  1 

OOt  3 

C‘ 

CONTINUE 

006  4 

NO=NO-  1 

0 0 6 5 

GO  TO  •;  1 0 , 1 2 . 1 4 , 1 6 . 1 8 . 2 0 ) , I N 0 U 

0 0 6 6 

10 

N07=NA 

0 06  7 

I FC  1 FLOG  . LT  . 2 1 )N07  = 0 

OOtS 

NB=NB+NA7 

0 0 6 9 

DO  11  J=1>JCHT 

0 0 ? 0 

X P 7 C J > ■-  X P C J > 

00?  1 

i 1 

NS7C  J ) = HSC  J ) 

0072 

GO  TO  6? 

0 07  3 

12 

N06  = NO 

0 0 7 4 

I FC  I U3 . LE  -6  >1 FL0G»0 

0 0 7 5 

IFC  IFLOG . LT  .21  )N06  = 0 

0 0 7 b 

N8=NB'*-N06 

0 0 ? 7 

DO  13  J= 1 , dCNT 

007  8 

XP6C  J > = XPC  J > 

0 0 7 9 

13 

N36C  J ) = HS<  J ) 

0080 

GO  TO  67 

0 08  1 

14 

N05=N0 

008  2 

IFC  IFLOG  . LT  21  )N05*0 

008  3 

NB=NB+N05 

0084 

DO  15  1 , J CHT 

0 0 8 5 

XP5C  J > = XPC  J > 

0 08  6 

15 

NS5C  0 > = NSC  ,1  > 

0 08  7 

GO  TO  67 

008  8 

16 

N04=N0 

008  9 

IFC  IFLOG  LT  .21  :)N04  = 0 

0 0 9 0 

NB=NB+N04 

009  1 

DO  1 7 = 1 , JCNT 

0 09  2 

XP4''  0 > = XP<  J ) 

009  3 

17 

NS4C  ..1  > = NS<  J ) 

0094 

GO  TO  67 

0 0 9 5 

18 

H03=N0 

009  6 

IFC  IUi3  LE  .-3>IFL0G«0 

0097 

IFC  IFLOG  . LT  .21  .>N03  = 0 

009  8 

NB  = NB-»-N03 

0099 

DO  19  d=l , JCHT 

0 100 

XP3C  J ) = XPC  J > 

0 10  1 

19 

MS3C  J > = HS<  J ) 

0 ! 02 

GO  TO  67 

0103 

20 

N02=N0 

0 104 

IFC  IFLOG  LT  .21  )H02*0 

0105 

NB=N6+N02 

0 106 

DO  21  J»1 idCNT 

0107 

XP2C  J ) = XP<  J ) 

0 108 
0109 

21 

67 

NS2<  J >*HS<  J ) 

CONTINUE 

0110 

IW2«2 

0111 

68 

CONTINUE 

0112 

0113 

0114 

0115 

IFC  N07  LE  .0  )C0  TO  25 

L-0 

DO  24  J»1,N07  ..i< 

KJ-NS7(J> 

0116 

DC  J >*0 . 0 

0117 

DO  22  K«1 >KJ 

0118 

L»L  + 1 
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0 119 

C'<  J >-[>(  J >+S(  L > 

0 1 2 0 

2 2 

COhT IHUE 

0 1 2 1 

8t>EC  = K J 

0 122 

D'!  J >-C;<  J >/6t>E  C 

0123 

B7(  J >=B7<.  u > + S(  22  )*[H  J ) 

0 12  4 

DO  2 3 I * 1 , J 

0 125' 

23 

A7(  I - J > = 87<  1 ..  J )'*-D<  1 :)*D<  d ) 

0 !2t 

24 

EOMT IHUE 

0 12  7 

'-.e 

CONTINUE 

0 12S 

1 F( H«6  . LE  . 0 JGO  TO  2 9 

012  9 

L=0 

0 13  0 

DO  28  J= 1 , NAB 

0 13  1 

D < J > = 00 

0 13  2 

Ko  = NS6>:  J > 

0 13  3 

DO  2 6 K = 1 - K J 

0 13  4 

L = L+  1 

0 135 

Di'  J > = D(  J -■  + s>;  L ; 

0 1 3(s 

26 

CONTINUE 

0137 

ADEC=Kd 

0 138 

D<  J > = D<  J )/ADEC 

0 139 

B6<  J >=e6<;  J >+s<  22  >*D':  j > 

0 140 

DO  27  I = 1 . J 

0 14  1 

27 

A6(  I ; J ) = A6(  I , J >*-D<  I )*D<  J > 

0 142 

28 

CONTINUE 

0143 

29 

CONTINUE 

0 144 

IF( NA5 . LE . 0 )G0  TO  33 

0 145 

L = 0 

0 14  6 

DO  32  J=1,NA5 

0 147 

i‘<.  J > = 0.0 

0 148 

K.J  = NS5';  J > 

0 149 

DO  30  K>1-KJ 

0 150 

L=L  + 1 

0 15  1 

D< J > = D<  J ) + S<L  ) 

0 15  2 

30 

CONTINUE 

0 153 

ADEC=Kd 

0 154 

DC  J > = D<:  J >/ADEC 

0 15  5 

B5(  J > = B5<:  J ) + S(  22  )*DC  J > 

0 15  6 

DO  3 1 1 = 1 , J 

0 15  7 

31 

A5(  I , J ) = A5<  I . J >4-D(  I >*D(  J > 

0 158 

32 

CONTINUE 

0 159 

33 

CONTINUE 

0 160 

IF( NA4 .LE  O >G0  TO  37 

0 16  1 

L=0 

0 162 

DO  36  J= 1 , NA4 

0 163 

DC  J > = 0 . 0 

0 164 

KJ  = NS4( J ) 

0165 

DO  34  K=1  >(CJ 

0166 

L=L  + 1 

0 167 

DC  J ) = DC  J >+SC  L > 

0 168 

34 

CONTINUE 

0 169 

ADEC=Kd 

0 170 

DC  J ) = DC  J :)/ADEC 

0171 

B4(  J >=B4C  J ) + SC  22  >*DC  J > 

0172 

DO  35  1= 1 , J 

0 173 

35 

A4C  I . J ) = A4C  I , J ) + D(  I )*D(  J ) 

0174 

36 

CONTINUE 

0175 

37 

CONTINUE 

0 176 

IFC  NA3 .LE .0  )G0  TO  41 

0177 

L=0 

0178 

DO  .40  J=  1 , NA3 

0 1 7'9 

O’'  J ) = 0 , 0 

0 16  0 

KJ  = H83'::  J > 

0 18  1 

DO  38  K=1 ,KJ 

0 18  2 

L=L+  1 

0 1 8 3 

0(J)~0>:J)+S(L) 

0 184 

38 

CONTINUE 

0 18  5 

8DEC=K J 

0 188 

D ( J > = D < J > / A D E C 

0 187 

B3(  J >=e3( J ) + S( 22  )*DC  J ) 

0 188 

DO  39  I = 1 , J 

0189 

39 

A3(  I . J >=A3(  I .•  J )+D(  I )*D<  J ) 

0 190 

40 

CONTINUE 

0 1 9 1 

41 

CONTINUE 

0 19  2 

I F(  HA2  . LE  . 0 )G0  TO  6 9 

0 19  3 

L = 0 

0 194 

DO  4 4 ..1=  1 , NA2 

0 195 

D<  J :>  = o . 0 

0 19  8 

K..l  = NS2<;j  > 

0 197 

DO  42  K=1 ,KJ 

0 198 

L=L  + 1 

0 199 

DC  J ) = D( J > + S<  L > 

0 20  0 

42 

CONTINUE 

0 201 

ADEC=KO 

0 20  2 

D C J > = D C J ) / A D E C 

0 20  3 

B2(  J ) = B2<  J ) + S(  22  )*D(:  J ) 

0204 

DO  4 3 I = 1 . J 

0 2 0 5 

43 

A2(  I .■  J )*A2(;  I .■  J .'»  + D(  I >*D(  J ) 

0 2 0 8 

44 

CONTINUE 

0207 

89 

CONTINUE 

0 20  8 

RETURN 

0 2 0 9 

END 

0 2 1 0 

END* 
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SULTMl  T=00'>04  IS  ON  CR00002  USING  OOOll  BLKS  R = OOSl 


0 00  1 
O00  2 
0 00  3 

■;>  />  n 4 
0005 

0 0 0 i 
000? 
0 o 0 6 
0 0 0 S 
0 0 1 0 
0 0 1 I 
n >;)  1 2 

0 0 i 3 
0 0 1 4 
0015 

001b 

0 017 
0 018 
001  S 
0 0 2 0 
002  1 
002  2 
002  3 
0 0 2 4 
0 02  5 
0 02  6 
002  7 
002S 
0 02  9 
0030 
0031 
0032 
0 03  3 
0034 
0 03  5 
0036 
003  7 
003  8 
0039 
0040 
004  1 
004  2 
0 04  3 
0044 
0 04  5 
0 04  6 
0 04  7 
0048 
004  9 
0050 
0 05  1 
0052 
005  3 
0054 
0055 
0056 
005  7 
005  8 


FTN4  L 

SUBROUTINE  WLTM3 

C ULTN3  IS  USEU  TO  COMPUTE  THE  FIT  FOR  UEIGHTING  FACTORS 

C FOR  DATA  FROM  MAG  TAPE. 

COMMON  F . S( 22  ) . I T( 6 ) 

COMMON  A7<4;4),B7(4).XP7<4)/NA7,NS7<4>,SB7,BA7<'4) 

COMMON  A6(  5 ..  5 ) , B6<  5 >/  XP6<  5 ).  NA6>  NS6<  5 ; , SB 6 . B A6 ( 5 j 

COMMON  A5':  5 , 5 ) , B5<  5 XP5(  5 NAS,  NS5(  5 ) , SB5.  BA5(  5 ) 

COMMON  A4<.6,6.),B4(6),XP4(6),Nh4,NS4(6),SB4,BA4';6> 

COMMON  A3<'8.6),B3(8),XP3(.8),HA3,NS3<8),SB3,BA3(8> 

COMMON  A2<  1 1 , 1 1 ),  62<  1 1 ) , XP2(  1 I >,  NA2,  NS2<  1 1 ) , SB2..  Bh2(  1 1 ■> 
COMMON  P , NB , MN , I S, MS, I ARRY 
COMMON  Ihil,  IU2,  1U3,  IU4,  IU5,  IlilZ 
COMMON  WT(21) 

COMMON  XC,YC,SN,SX,SY,SXX,SYX,DA,DS 
DIMENSION  A<  1 1 , 1 1 ) 

WRI  TEC  6-  99  > I T( 4 ) , I T( 3 ), P 

99  FORMaTC IHO,  "AT",  13,  “ ! 12,  " U I T H " , F 1 0 . 0 , " POINTS”) 

DO  69  L= 1 , I U5 

GO  TO  C 1,3, 5, 7, 9, 11), L 

1 CONTINUE 

IFCNA? .LE  .0  )G0  TO  69 
DO  2 J =1  ..  NA7 
BA7C  •)  )=B7C  J > 

DO  2 1=1 , J 
A<  I , ..I  ) = A7<  I , J ) 

2 AC  J , I ) = A7<  I , J ) 

NG=7 

MTS  = 4 

CALL  W3SUB(A,B7,XP7,NA7,NS7,SB7,BA7,NG,MTS) 

I FC  NG . EO  . 7 )G0  TO  69 

NB=NB-NA7 

NA7=-1 

GO  TO  69 

3 CONTINUE 

IFC  NA6  LE .0  )G0  TO  69 
DO  4 J = 1 , NA6 
BA6C  J ) = B6<  J ) 

DO  4 1=1,0 
AC  I , J )=A6':  I , J ) 

4 AC  0 , I )=A6C  I , J ) 

NG=6 

MTS  = 5 

CALL  lil3SUBC  A,  B6,  XP6  , NAB  , NS6,  SB6  , BAS,  NG  , MTS  > 

IFC NG. EQ . 6)G0  TO  69 
N6=NB-NA6 


NA6=-1 
GO  TO  69 

5 CONTINUE 

IFC  NAS  LE . 0 )G0  TO  69 
DO  6 0=1, NAS 
BA5( 0 ) = B5C  J ) 

DO  6 1=1,0 
AC  I , 0 )*A5C  1,0) 

6 ACO , I )=A5<  I ,0  ) 

NG  = 5 
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005  9 
0 06  0 
006  I 
0062 
0063 
0 06  4 
0 0 6 5 
0066 
006  ? 
006  8 
0 06  9 
0 0 7 0 
007  1 
0 07  2 
0 07  3 
0074 
0075 
0076 
0077 
00  73 
007  9 
008  0 
008  1 
008  2 
0 08  3 
008  4 
0 0 8 5 
0 08  6 
008  7 
0088 
0 08  9 
009  0 
0 09  1 
0 09  2 
0 0 9 3 
0 09  4 
0 09  5 
0096 
009  7 
0098 
0099 
0 1 0 0 
0 10  1 
0 102 
0 103 
0104 
0 105 
0 106 
0 107 
0 108 
0 109 
0 110 
0 111 
0112 
0113 
0114 
0115 
0116 
0117 
0119 


MTS  = 5 

CALL  W3SUB<  A..  85,  KP5  , NA5  , HS5,  SB5  , BA5,  NG  , MTS  > 

IF( NG  £0  5 )G0  TO  69 

NB=HB-NA5 

NA5=-l 

GO  TO  69 

7 CONTINUE 

I F( NA4  LE  . 0 )G0  TO  69 
DO  8 J =1  • NA4 
B A 4 < 0 ■>  = B 4 < .1  > 

DO  8 1=1.0 

I , j > = fl4t'  I , 0 > 

8 A<  J , I ;)  = A4«'  I , J > 

NG=4 

MTS  = 6 

CALL  til3SUB<  A,  64,  XP4  , NA4  , NS4,  SB4  , BA4,  NG  , MTS  ) 

I F( NG . EG  4 >C0  TO  69 

NB=NB-NA4 

NA4=-1 

GO  TO  69 

9 CONTINUE 

IF(Nih3  LE.OIGO  TO  69 
DO  10  0=1, NA3 
B A 3 < J > = B 3 <•'  J > 

DO  10  1=1.0 
A < I , J ■)  = A 3 ( I • 0 > 

10  A<  0 , n = A3<:  I , 0 > 

NG=3 

MTS  = 8 

CALL  W3SUB<  A, 6 3, XP3 , NA3 , NS3, S83, BA3, NG , MTS  ) 

I F(  HG  EG  . 3 .'GO  TO  69 

HB=HB-NA3 

NA3=-1 

GO  TO  69 

11  CONTINUE 

IF<  NA2  LE  O >G0  TO  69 
DO  12  0=1 ,NA2 
BA2<  J )=B2<  J ) 

DO  12  1=1,0 
A< I , 0 )=A2<  1,0) 

12  A( J , I ) = A2C  1,0) 

NG=2 

MTS=  1 1 

CALL  U3SUB<  A,  62, XP2,NA2, NS2, SB2,BA2, NG, MTS) 

IF( NG . EG . 2 )G0  TO  69 

Ne=NB-NA2 

NA2=-1 

69  CONTINUE 

RETURN 
END 

SUBROUTINE  U3SU6( A , B, XP , NA , NS , SB, BA, NG , MTS  ) 
DIMENSION  A(  1 1 , 1 1 ),  B<  MTS  ),  XP<  MTS),  NS<nTS),BA<:MTS  ) 
M = NA 
SB  = 0 , 0 


A11»A<  1,  1 ) 

IF( All  EQ  .0  )G0  TO  68 
DO  1 1=2, M 

1 A<  1 , I ) = A<  I,  1 )/Al  1 

6A(  1 )>BA<  D/All 


yr 
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DO  5 J=2,M 

J 1 = J - 1 

DO  3 I=J.M 

A S = 0 . 0 

DO  2 K = 1 . J 1 

AS  = HS+ A(  I . K K . J > 

A<  I ^ J > = A<  I J )- AS 

IF(  1 -GT.J  I)  — A(  ) 

CONTINUE 

BS  = 0 . 0 

DO  4 K = 1 , J 1 

BS  = eS  + A(  J ,K  >*BA( K) 

A j J = A ( .J  ^ J ') 

I F(  AJ J . EQ . 0 )GO  TO  68 
eA(  J > = < BA<  J >-BS  >/A  J J 
N 1 = fi  - 1 
DO  7 I = 1 , M 1 
BS=0 . 0 
M I = M - I 
M 1 1 =MI +1 
DO  6 J=MI  1.  M 
BS=BS+A(  MI. J )*BA< J > 

BA(  MI > = 6A(MI )-BS 

CONTINUE 

WRI  TE(  6.  -JS  /NG 

FORMATS  IHO,.  “FOR  GROUPS  OF". 13) 
WRITE<6.93)<NS<J).J  = 1,M:) 

FORMATt-  IH  . ,1116) 

WRI  TE<  6.  9 7)<  XP<  J ).  J = 1 . « ) 

FORMAT< IH  . "X"  . 1 1F6 . 1 ) 

WRITE(6.96>(BA<J  ).J  = 1.M) 

FORMAT*;  IH  . "Y“  . 1 1F6  . 3 ) 

DO  8 I =1 . M 
SB=SB+BA(  I ) 

CONTINUE 
MRI  TE(  6. 85  )30 

F0RMAT<1H  ."SUM  OF  WEIGHTS  =",F8.5) 

GO  TO  68 
CONTINUE 
URI  TE<  6.  94  )NG 

FORMATS  IHO.  "FOR  GROUPS  OF". 1-3."  MATRIX 

NG=-  1 

CONTINUE 

RETURN 

END 

END$ 
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JSHIS  PAUi  1=  ‘ _ 

„ . jjj  C 

fROU  Gux 1 ■ 


fRACTlftlBI 


$WLTM4  T=00004  IS  OH  CR00002  USING  00012  BLKS  R=00?7 


OC'O  1 
0 0 0 2 
0 0 0 3 
0 O 0 4 
0 0 0 5 
0 0 0 1. 
0 0 0 7 
0 0 0 3 
0 0 0 9 
0 0 1 0 
OOl  1 
001  2 
0 0 1 3 
0 0 1 4 
0 0 1 5 
0 013 
0017 
0 0 1 S 
0019 
002  0 
0 0 2 1 
002  2 
0 02  3 
0 02  4 
0 02  5 
0 02  6 
0 02  7 
002  8 
002  9 
0 0 3 0 
■j  0 3 1 
003  2 
0 03  3 
0034 
003  5 
0 0 3 C 
0 0 3.7 
0 03  3 
0 03  9 
0040 
004  1 
0042 
0 04  3 
004  4 
0 04  5 
004  6 
0047 
0043 
0 04  9 
0050 
0 051 
005  2 
005  3 
0054 
0055 
0056 
0057 
005  8 
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PROGRriM  ULTN4.3 

c + + 

C UILTM4  IS  H PROGRAM  WHICH  NORMftLIZES  THE  RESULTS  OF  ULTM3, 

C PLOTS,  REPORTS,  HND  PUNCHES  A TAPE  OF  THESE  RESULTS. 

C IT  CAN  ALSO  RESCHEDULE  WLTMG  FOR  A REPEAT  WITH  A DIFFERENT 

C SLIDE  FACTOR. 

j;  * * * * :*;f  * ♦;*  .d*  .■*  * ******♦*#**♦♦*♦♦*♦  .■*♦;♦♦♦+**♦♦  ;(i  ♦♦♦  ♦**:(!  .+ 

COMMON  F , S( 22  >, I T( 6 > 

COMMON  A7<:4,4>,e7(4>,XP7<4>..NA7,NS7(4;'.SB7,BA7(4:' 

COMMON  A6<.  5 , 5 > , B6<  5 .1,  XP6<  5 >,  NA6  , NS6(  5 > , SB6,  Bm6<  5 ) 

COMMON  A5!.  5 , 5 > , B5(  5 >,  XP5<  5 .).•  NA5  , NSSI  5 ',' , SB5,  BA5(  5 ) 

COMMON  A4<  6 , 6 ) , B4<  6 ),  XP4<  6 NA4  , NS4<  6 >,  SB4,  BA4C  6 - 

COMMON  A3<  8 , 8 > , B3(  8 ),  XP3<  8 ),  NAS  , NS3<  8 ) , SB3..  BA3<  8 ) 

COMMON  A2';  1 1 , 1 1 > ,e2<  1 1 ) , XP2<  1 1 > , H A 2 , NS  2 < 1 1 > , SB  2 , B A2  < 1 1 > 
COMMON  P, NB,MN, IS, MS, lARRY 
COMMON  I M 1 , IW2 , I W3 , IW4, I W5 , I W2 
COMMON  WTi:  2 1 > 

COMMON  XC,YC,SN,SX,SY,SXX,SYX,DA,DS 
DIMENSION  NMG<  3) 

DATA  NMG/2HWL , 2HTM , 2HG  / 

CALL  EXEC<  3 , 1 004B  ) 

WRI TE( 4, 99  >NB 
99  FORMAT-:  I 3.1 

DO  7 M = 1 , I W 5 

GO  TO  < 1 , 2, 3, 4 , 5 , 6 ) , H 

1 I F(  NA7  LE  . 0 :>G0  TO  7 
NG=7 

MTS  = 4 

CALL  W4SUB<XP7,NA7,NS7,SB7,BA7,NG,MTS) 

GO  TO  7 

2 I F(  HA6  . LE  . 0 IGO  TO  7 
NG  = 6 

MTS  = 5 

CALL  W4SUB(  XP6,NA6, NS6, SB6,BA6, NG, MTS) 

GO  TO  7 

3 IF( NA5 .LE .0  )G0  TO  7 
NG=5 

MTS  = 5 

CALL  U4SUB<  XP5,NA5, NS5, SB5,BA5, NG, MTS) 

GO  TO  7 

4 IF(  NA4  . LE  . 0 )G0  TO  7 
NG=4 

MTS  = 6 

CALL  W4SUB<  XP4 ,NA4 , NS4, SB4 , BA4, NG , MTS  ) 

GO  TO  7 

5 I F( NAS . LE . 0 )G0  TO  7 
NG  = 3 


6 


7 


MTS  = 8 

CALL  «4SUB<  XP3,NA3, NS3, SB3,BA3, NG, MTS) 
GO  TO  7 

IF(NA2.LE.0)G0  TO  7 

NG>2 

MTS=11 

CALL  U4SUB( XP2,NA2,NS2, S62,6A2, NG,MTS> 
CONTINUE 

CALL  EXEC< 1 1, I T) 
WRITE<6,98)IT<4),IT<3),IT<2) 


53 


F OR  MAT  C 1 HO..  “ANALYSIS  CO  M PL  E T E D " , I 3 . 


CALL  EXEC<3.100 
I F(  I 04  . EQ . 1 >G0 
I F<  I 03  . GE  0 ,'GO 
103=103+1 
I F(  I 03  EQ  . 0 )I  (.14 

n.i  1 = 0 
102=1 

ZERO  OUT  ARRAYS 
00  8 1=1.11 
DO  8 J = 1 . 1 1 
A2<  I . J :'  = 0 . 0 


3 . 1 0 0 4 B ) 

1 >G0  TO  88 
0 G 0 TO  8 8 

, 0 >104=  1 


0LTM2 


GT  . 4 .-JGO 


A2<  I . J :'  = 0 . 0 
B2(  I ; = 0 . 0 
NS2<  I >=0 

IF(  I GT  . 8 . OR . J . GT . 8 )G0  TO  8 
A3(  I , J > = 0 . 0 
B 3 ( I ) = 0 . 0 
N S 3 ( r;  = 0 

I F(  I . GT  . 8 . OR  . J . GT  . 8 :)G0  TO  8 
A4(  I , J > = 0 . 0 
B4(  I .■>  = 0 . 0 
NS4<:  I >=0 

IF(  I GT  . 5 OR  . J . GT  . 5 .)G0  TO  8 

A 5 ( I , J > = 0 0 

B5<  I ) = 0 . 0 

HS5(  I .:'=0 

A8<:  I • J ) = 0 . 0 

B8(  I > = 0 . 0 

NS6<  I )=0 

I F(  I . GT  . 4 . OF:  . J GT  . 4 .-JGO  TO  8 
A7(.  I , J > = 0 0 
B7t:  I :'=o . 0 
NS7<;  I > = 0 
8 CONTINUE 

P = 0 . 0 

CALL  EXEC< 1 0. NMG  ) 

68  CONTINUE 

STOP 
ENO 

SUBROUTINE  04SUe<  XP  . NA. NS, SB . BA . NG . MTS  > 

DIMENSION  XP<  MTS  ),  NS<  MTS  BA(  MTS) 

CALL  PLTLU(  10  ) 

CALL  SFACK  15,10.) 

CALL  LLEFT 

CALL  PLOT<  0 .0,0.0,-!) 

CALL  PLOT<  3 5,1.0,3) 

S6  = ABS<  SB  ) 

DO  1 1=1, NA 

BA(I)  = BA(:i>/(SB*FLOAT<NS<I))) 

XV  = XP(  I )/2  +3.5 
YP  = BA<  I )*50  +1  . 

IF( YP. LE . 0. 2)YP=0. 2 
I F<  YP . GE . 8 8 )YP-8 . 8 
CALL  SYMB< XV, YP, 0 . 1 4, NC, 0 0, -1  ) 

1 CONTINUE 

IrfRI TE( 8, 88>NG 

88  FORMAT( IHO,  "NORMALIZED  WEIGHTS  FOR  CROUPS  OF",  13) 

URITE<6,88)<NS(I),I  = 1,HA) 

88  FORMAT< IN  ,"#",1118) 


1 4, NC, 0 0, -1  ) 


CROUPS  OF",  13  ) 


XUIS  PAGE  , 

FRUId  Uti  if' 


^.ALITY  PlUaiUkiJM 
i,:,U)DDC  


0 1 1 •? 
0 12  0 
0 12  1 
0 1 2 2 
0 123 
0 12  4 
0 125 

0 12b 

0 127 
0 1 2 8 


WRI  TE(  lew  ?7  I ),  1 = 1 , NA  ) 

■?7  FORMAT-:  IH  . “K"  , 1 1F6  . 1 > 

WRI  TE(  6.  96  eA<  I I =1  / NA  ) 

96  FORMAT-;  IH  , " Y " , 1 1 F 6 , 3 ) 

yRI  TE<  4.  95  >C  NG  . NSC  I ) / XP-:  I >.  BA<  I >,.  I = 
95  F OR  M AT  C I 2 . " , " , I 3 , % F5  . 1 . - F 9 . 5 > 

CALL  LLEFT 
RETURN 
END 
END« 


.0^- 


1 / NA  ) 
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SiiiLTM?  T=00o()4  IS  ON  CRC'0002  USING  000  IS  BLKS  R=0142 


'"'0  1 
0 0 2 
0 00  3 
0 0 0 4 
0005 
OoOt 
‘0  0 0 7 
0 0 0 S 
0 0 0 9 
0 0 1 0 
0 0 1 1 
0 0 i 2 
0013 
0 0 1 4 
0015 
0 0 1 6 
0 0 1 7 
0 0 1 S 
0 0 1 9 
0 0 2 0 
0 02  1 
0022 
0 02  3 
002  4 
0 02  5 
0 02  6 
0 02  7 
002  8 
0 02  9 
0 0 3 0 
0 031 
003  2 
0033 
003  4 
003  5 
O'  o 3 6 
0 03  7 
003  8 
003  9 
004  0 
004  1 
0042 
004  3 
004  4 
004  5 
0046 
0 04  7 
0 04  8 
004  9 
0050 
005  1 
0 05  2 
0053 
005  4 
005  5 
0056 
0 05  7 
0058 
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SUBROUTINE  ULTM5 

C ULTMS  IS  H PROGRAM  UHICH  WILL  REPORT  THE  NECESSARY 

C INFORMATION  ON  PLOTS  AND  DRAW  THE  GRID  IF  RE&UIRED 

C 

COMMON  F . S<  22  I T( 6 ) 

COMMON  A7(4,4),e7<4>,XP7(4;,NA7,NS7<4),SB7.BA7<:4> 

COMMON  A6^5,5>,e6(5),XP6(5).NA6<HS6(5),SB6.BA6<5) 

COMMON  AS*:  5 . 5 :>  ,B5<  5 XP5(  5 NA5  < NS5(  5 > , SB5-  BA5',  5 ) 

COMMON  A4<6.6>,B4(6>,XP4<6).NA4.NS4(6>,SB4.Ba4>;6> 

COMMON  A3(  8 - 8 :>  - B3<  8 ><  XP3(  8 )i  NA3  - NS3<  8 ) , SB3.  BA3(  8 ) 

COMMON  A2(  1 I , 1 1 ),  B2(  1 1 ) < XP2(  1 1 > , NA2.  NSBC  1 1 > , SB2..  Ba2<  1 1 :> 
COMMON  P>NB.MN,IS.MS. lARRY 
COMMON  I M 1 , IW2 . I W3 , IU4, I W5, I U2 
COMMON  WT( 2 1 > 

COMMON  XC , YC, SN, SX, SY, SXX, SYX . DA. DS 
DIMENSION  NA<2>.NS(2>,NPi:2:i.NFS<2).NFL<2> 

DIMENSION  MON<  2>.NCH( 2). NSEC<  2> 

DATA  NA/2HAV. 2HG=7 . NS/2HSD . 2HV=/ 

DATA  NP,''2H#P,  2HTS/ 

DATA  HFL/2HFI , 2HLE/ , NFS/2HFS . 2H=  / 

DATA  NSEC/BHSE . 2HC  / 

DATA  NCH/2HCH.2Hi  / 

NY=54475B 
NX=54053B 
CALL  PLTLU(  10  ) 

CALL  SFACT( 15.10,  > 

CALL  LLEFT 

CALL  PLOT( 0 . 0 . 0 , 0. - 1 ) 

I F(  I W1  EQ  . 0 )G0  TO  5 
CALL  PLOT( 0 . 5. 5 . 0. 3 ) 

CALL  PL0T<9.5.5. 0.2) 

CALL  PLOT<  5 . 0. 0 . 5. 3 ) 

CALL  PLOT< 5 . 0 . 9 5. 2 ) 

CALL  DASHCO .5 . 0 . 5. 0 .5.0 . 5. -1  ) 

CALL  DASH< 0 5. 0 . 5. 9 . 5.9 . 5.  1 ) 

TO  69 


CALL  DASHC  0 5.0  5.0 

CALL  DASH<  0 5.0,5.9 

IF<  IWl  . LE  . 1 )G0  TO  6 

D=SX*SX-SN*3XX 

A=( SX*SY-SN*SYX )7D 

B = ( SX’t'SYX-SY*SXX  )/D 

X=-4 . 0 

CONTINUE 

Y=A*X-*  6 

Z = ABS<  Y) 

IF<  Z ,LE  . 4 .75)G0  TO 

X=X+  5 

GO  TO  1 

CONTINUE 

X=X+5 . 

Y=Y+5 . 

CALL  PL0T<X.Y,3) 

X = 4 . 0 
CONTINUE 
Y*A*X+B 
2=ABS<  Y) 

IF<  Z LE  4 75>G0  TO 
X=X-0 . 5 


56 


GO  TO  3 
COMT  INUE 
X=X+5 . 

V=Y+5 . 

CALL  FL0T<X.Y,2> 

DATT=DAXSN 

DSTT=SQRT<<DS/SN)-DATT*DATT) 

CONTINUE 

CALL  LLEFT 

CALL  PLOT<  0 . , 0 . . -1  ) 

IDAY  = IT<  5 > 

I YEAR= ITC  6 ) 

CALL  DATE< IDAY,MOH,  lYEAR  ) 

DAY= IDAY 
YEAR=I YEAR 
FILE=IWZ 
CHN=  IS 
SEC=MS 

IF(  MS. EQ . 0>SEC  = 0 . 5 
FS=5 . /<  XC*30 . ) 

CALL  NUMB<  1 . 0 . 9 . 0.  0 . 1 4,  DAY  , 0 . 0 -1  > 

CALL  SYMe<1.56.9,0,0.14..  MON.0.0,3) 

CALL  NUNe( 2 , 25 , 9 . 0 , 0 . 14. YEAR / 0 . 0. -1  > 
CALL  SYMB<  1 . 0 . 8 . 5.  0 . 1 4.  NCH  . 0 . 0 ..  3 ) 

CALL  NIJMB(  1 . 58 . 8 . 5. 0 . 14  . CNN..  0 . 0 , - 1 ) 
CALL  NUMe<  2 . 25 , 8 . 5 . 0 , 14  . SEC..  0 , 0 . 1 ) 

CALL  SYMB'C  2 . 85. 8 . 5 . 0 . 14  . NSEC  . 0 . 0.  3 ) 
CALL  SYMB<3.40.8.5,0. 14.NFL.0.0.4> 

CALL  NUMB<4.10.8,5.0. 14. FILE. 0,0. -1) 

I F<  I W1  . EQ  . 0 >G0  TO  69 

CALL  SYMB<1 .0.8. 0.0. 14. NFS. 0.0. 3) 

CALL  NUMB <999. 0.999. 0.0.  14.FS.0.0.  1) 
CALL  SYMB<1 .0.7.5.0.14.NY.0.0.2> 

CALL  NUMB<  999 . 0.  999 . 0 . 0 . 14  . A.  0 . 0.  3 
CALL  SYMB<  999 . 0. 999 . 0 . 0 . 14 . NX . 0 . 0 . 2 ) 
CALL  NUMB<  999 . 0.  999 . 0 . 0 . 14 . B. 0 . 0. 3 ) 
CALL  SYMB< 1 , 0 . 7 . 0. 0 . 14. NA. 0 , 0 . 4 ) 

CALL  NUMB<999.0.999.0.0.14.DATT.0.0.3) 
CALL  SYMB<3.0.7.0.0.14.NS.0.0.4) 

CALL  NUMB<  999 . 0.  999 . 0. 0 . 1 4 . D S TT . 0 . 0 . 3 ) 
CALL  SYMB< 1 . 0 , 6 . 5. 0 . 1 4. NP. 0 . 0 . 4 ) 

CALL  NUMB<1.60.6.5.0.14.SN.0.0.-1) 

CALL  LLEFT 
CONTINUE 

IF<  lUl  EQ  . 1 >1W1:==2 
RETURN 
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